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ABSTRACT 

We mapped the MSX dark cloud G084.81-01.09 in the NH3 (1,1) - (4,4) hnes 
and in the J = 1-0 transitions of ^^CO, ^^CO, C^^O and HCO+ in order to study 
the physical properties of infrared dark clouds, and to better understand the 
initial conditions for massive star formation. Six ammonia cores are identified 
with masses ranging from 60 to 250 Mq, a kinetic temperature of 12 K, and a 
molecular hydrogen number density r;,(H2) ~ 10^ cm~^. In our high mass cores, 
the ammonia line width of 1 km s~^ is larger than those found in lower mass cores 
but narrower than the more evolved massive ones. We detected self-reversed 
profiles in HCO^ across the northern part of our cloud and velocity gradients 
in different molecules. These indicate an expanding motion in the outer layer 
and more complex motions of the clumps more inside our cloud. We also discuss 
the millimeter wave continuum from the dust. These properties indicate that 
our cloud is a potential site of massive star formation but is still in a very early 
evolutionary stage. 

Subject headings: stars: formation - ISM: molecules - ISM: kinematics and dy- 
namics - ISM: individual (G084.81-01.09) 

1. INTRODUCTION 

Infrared dark clouds (IRDCs) are identified as small areas with high extinction against 
the bright diffuse mid-infrared background along the Galactic Plane. These "flecks" in the 
sky were first discovered by the 15 /xm ISOGAL survey (the inner Galactic disk survey of 



^Purple Mountain Observatory, Chinese Academy of Sciences, Nanjing 210008, China; 
shbzhang@pmo .ac.cn 

^Institute for Astronomy, 2680 Woodlawn Drive, Honohilu, HI 96822, USA 

•^Max-Planck-Institut filr Radioastronomie, Auf dem Hitgel 69, 53121 Bonn, Germany 



-2- 



the Infrared Space Observatory, iPerault et al.l Il996l ) and later cataloge d and studied us - 
ing the 8.3 jUm survey data of the Mid-course Space Experiment (MSX, lEgan et al.lll998l ). 
Egan et al.l found about 2000 IRDCs and note that these objects are cold and dense molecu- 



lar cores. These resu lts were refined by later molecular li ne observations (ICarey et al.lll998 



Teyssier et al. 2002) and radio continuum observations ( Carey et al. 200ol ). These authors 
concluded that these clouds have low kinetic temperatures of 10-20 K and a gas density over 
10^ cm~^. Their small line widths of 1 - 3 km s~^ can be related to a very early phase of star 
formation. Therefore, the IRDCs are good candidates for pre-protostellar cores and may be 
the primary sites for future massive star formation, which plays an important role in the 
evolution of our galaxy. 

However, IRDCs remain mysterious in at least some aspects. A controversial issue 
being discussed is whether the cloud evolution is quasi-static or whether it is undergoing 
a dynamical process caused by lack of equilibrium ( iBergin fc Tafallall2007l ). Understanding 
the role of supersonic turbulence in dark clouds is important in settling this argument. 
Evidence from current observations of IRDCs seems to support both v iews. Most IRDC 
cores are supported by turbulent pressure but appear to be virialised (jPillai et al.l l2006l ) . 
Furthermore, the origin of the turbulence in the molecular clouds and the transition from 
the turbulent cloud to the quiescent core regime is still unclear. Also, more systematic studies 
need to be done to compare them with other relatively nearby massive star-forming clouds. 
The answers will provide a critical element in our general understanding of star formation. 

In this study, we have obtained spectroscopic data in several mo lecules towards the 
dark cloud G048. 81-01. 09, a source in the catalog of ISimon et al.l ( l2006l ). The cloud has an 
elongated morphology (N- S) with an extent of ~ 15 arcmin, and is situated i n a region o f 
high extinction {A^ > 20, ICamberesv et al.ll2002h of the dark cloud LDN935 JLvndslll962[ ) 
which separate s the W80 HII region into th e North America (NGC 700 0) and Pelican (IC 
5070) nebulae (IComeron fc Pasqualill2005l ). IComeron fc Pasqualil (120051 ) proposed that the 
ionizing source of the complex is 2MASS J205551. 25+435224. 6, an 05V star located behind 
LDN935. The IRDC is almost starless with no associated IRAS or MSX point sources. 
Previous molecular line observations indicated the presence of large amounts of gas with 
different vel ocity components and s igns of ongoing star formation in LDN935 around our 
dark cloud (iFeldt &: Wendkerl Il993l ). A single pointed CO observation towards our cloud 
conducted by Milman ( 19751 ) suggests however that the gas is not appreciably heated. 



There is no direct distance measurement for G084. 81— 01.09, but it is associated with 
dark cloud LDN935 for which there are a number of distance estimates. A radio continuum 



study carried out by lWendker et al.l (119831 ) gave a distance estimate of 500 pc. IStraizys et al. 



(119931 ) obtained a consistent result of 580 pc on the basis of photometric results of 564 
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stars toward three are as near the North Ame rica and Pehcan Nebulae complexes. New 
measurements done by ICersosimo et al.l (120071 ) derived a distance of about 0.7 kpc, which 
was obtained from radio recombination lin e observations at a frequency around 1.4 GHz and 
by applying a fiat rotation-curve model. IComeron fc Pasqualil ( 120051 ) obtained a distance 
of 610 pc to the ionizing star, 2MASS J205551. 25+435224. 6. Here we adopt a distance of 
580 pc in our calculation of physical parameters. 



2. OBSERVATIONS AND DATA REDUCTION 

2.1. Effelsberg 100 m Observations 

We mapped the infrared dark cloud G084.81-01.09 in the NH3 (1,1), (2,2), (3,3) and 
(4,4) transitions using the Effelsberg 100 m telescope of the Max-Planck-Institut fiir Ra- 
dioastronomie (MPIfR) during December 12 - 14, 2007 (see Table [1] for adopted line fre- 
quencies). The spectrometer backend employed was the 8192 channel AK90 auto-correlator 
which consists of 8 individual 1024 channel correlators of 20 MHz bandwidth each. Addi- 
tional data of NH3 (1,1), (2,2) and (3,3) were obtained on February 21 - 24, 2008 using the 
Fast Fourier Transform Spectrometer (FFTS), which has 16382 channels covering a band- 
width of 500 MHz. The data were taken in using frequency switching with a frequency throw 
of 7.5 MHz. The map covers the area with notable extinction in the MSX 8.3 //m image 
with a grid spacing of 30". The pointing was checked roughly every hour by observations of 
nearby continuum sources and was found to be accurate to better than 10". 

Continuum scans were used for calibrating the absolute flux density. The raw data were 
scaled from arbitrary units to main beam brightness temperature by applying gain-elevation 
corrections and the flux scale was se t using NGC 702 7 by assuming it to have a flux density 



of 5.4 Jy at our observing frequency lOtt et al.l (Il994j ). Our absolute calibration uncertainty 
is ±10%. 



2.2. PMODLH 13.7 m Observations 

The i^CQ (1-0), ^'^CO (1-0), C^^O (1-0) and HCO+ (1-0) observations of the 
MSXDC G084. 81— 01.09 (line frequencies shown in Tabled]) were carried out at the Purple 
Mountain Observatory Delingha (PMODLH) 13.7 m telescope from March to June in 2008. 
The 3 mm SIS receiver was used in double sideband mode. Using three 1024 channel AOS 
spectrometers with bandwidths of 145.43 MHz, 42.87 MHz and 43.3 MHz, the three CO lines 
were observed simultaneously. The source was mapped using position switched observations. 
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with the standard chopper wheel method for cahbration (jPenzias &: Burruall973l ). Standard 
sources were checked roughly every 2 hours. 

To compare with the NH3 data, our CO map covered the area mapped with NH3 and 
extended a few arcmins to the northeast. The grid spacing was 30" in the center and 60" 
in other regions. The HCO"'" map covered a region similar to that of the NH3 map with a 
typical grid spacing of 30". 

A summary of all observation parameters is provided in Table [H 



2.3. Data Reduction 



We used the CLASS software packag^ for the spectral data reduction. After discarding 
the bad scans, the spectra of each position were averaged and a polynomial baseline of orders 
1-3 was subtracted. Bad channels were excluded from the baseline fitting. 

The method "NH3(1,1)" in CLASS JBuIssou et allboosh was used to fit the hyperfine 



structure of the NH3 (1,1) spectra and to derive optical depths and line widths. The (2,2) line 
is so weak that only the main line could be dectected towards a few positions. Therefore, 
we simply fitted the (2,2) line with a single Gaussian to derive its main beam brightness 
temperature and an optical depth was not calculated in this case. Table [2] summarizes the 
fit parameters toward the cores detected in the data. No NH3 (3,3) or (4,4) lines were 
detected in our observations. 

The excitation temperature, rotational temperature, kinetic temperature and ammonia 
column density toward the cores are listed in Table [3l These physical parameters were 



derived using the standard formulae for NH3 spectra (IHo &: Towneslll983l ). 



The CO and HCO+ main beam brightness temperatures were derived from the antenna 
temperature, T^, using a main beam efficiency, r/mb? of 61%. A first order baseline was 
applied for the CO spectra and lst-3rd order baselines were used for the HCO^ spectra. A 
single Gaussian fit to the spectra provided the line velocity with respect to the local standard 
of rest (LSR) and the line width for later discussion. The fit results are given in Table |H 



^ CLASS is part of the GILDAS software suite available from 



http : //www . iram . f r/IRAMFR/GILDAS 
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3. RESULTS 

3.1. Identification and Morphology of the Cloud 

Fig. [T] shows the NH3 integrated intensity maps created from the velocity range -2.5 
to 4.5 km s~^ (which excludes the satellite lines). Our images of ammonia emission toward 
G084. 8 1—01.09 reveal extended, filamentary molecular emission that closely matches the 
morphology of the Spitzer mid-infrared extinction. Also, the ammonia maps (both NH3 (1,1) 
and (2,2)) allow us to identify two main molecular condensations: one in the northeast, 
containing ammonia cores 1, 2 and 6, and one in the south, consist ing of cores 3, 4 and 5 . 



All cores defined in the NH3 (1,1) map follow the standards given by lTachihara et al.l (120001 ) 
and the name of the cores were designated in descending order of their integrated intensities, 
as tabulated in Tabled The NH3 (1,1) core sizes range from 0.25 to 0.42 pc. 

The NH3 molecular emission also matches the 1.1 mm dust emission of Bolocam Galac- 



tic Plane Survey (BGPS. lAguirre et al.ll2010l ) remarkably well as shown in Fig. [2^. Further- 



more, with the blue cir cles marking the clump peaks of continuum emission identified by 



Rosolowsky et al.l (120091 ). we note that the ammonia cores are coincident with the peaks in 



1.1 mm emission. 

The integrated intensity maps of ^'^CO (1—0) and C^'^O (1—0) are shown in Fig. [2)d 
and [2t. The ^^CO and C^^O emission lines show an extended feature along north-south 
direction, with a sharp cutoff toward the northwest. The extent of C^^O emission is larger 
than that of both the ammonia and continuum emission. The ^^CO line, which tends to 
trace a lower gas density than C^^O, appears to be more extended. The strongest part of 
^■^CO emission coincides with the northern condensation (peak at the position of ammonia 
core 1) and shows a tail to the east, while the ^^CO emission in the southern part is fiat with 
no peak associated with that in the C^^O map. 

Fig. [2]i shows the integrated intensity of HCO"'" (1—0) overlaid on dust emission. The 
HCO"'" emission shows a different distribution compared to CO. Since the line is optically 
thick (see §3.31) . it is likely to trace the envelope surrounding a dense core, while the C^^O 
line is expected to trace the dense core itself. 

A candidate mas sive young stellar object (MYSO), G084.7847-01.1709, identified by 



Urquhart et al.l (120091 ) in their 6 cm VLA survey is marked with cross symbol on Fig. [TJ It is 
spatially coincident with the north condensation of our cloud. However, no compact source 
can be found near this object in the Spitzer 24 or 70 fim data. The lack of infrared emission 
could be due to the source being an extragalactic object rather than a Galactic MYSO. 
Further observations at centimeter wavelengths are needed to determine the spectral index 
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of emission and the nature of this object. 



3.2. Line intensities and kinetic temperature 

The excitation temperature, Tgx, of the NH3 (1,1) transition (Column Tex in Table [3]) 



was obtained from the optical depth via the relation 
T^b = ^[J(Tex) - J(Tbg)](l - e-^) 



J(T) 



fe^ - 11 



where Tmb and r represent the temperature and the optical depth from CLASS fitting pro- 
cedures and Tbg equals 2.7 K. The rotational temperature between the NH3 (1,1) and (2,2) 



by 


Ho & Townes ( 


1983) 


Tafalla et al. 


( 


2004 


)• 



9831 ). We then estimate the kinetic temperature using the expression of 



The typical kinetic temperature in the cores is about 12 K as seen in Table [31 No 
obvious difference in temperature appears among the si x cores. The kinetic temperature 
in the envelope is about 2 K higher than in the cores. iDesert et al. J2008h carried out a 
low resolution (~12 arcmin) survey at four (sub) millimeter wavelengths and derived a dust 
temperature of 8.5 K with a fixed emissivity law exponent of 2 toward our cloud, which is 
slightly lower than our gas temperature, probably because of the different beam size. 

The physical parameters of CO are calculated under the assumption of Local Thermo- 
dynamic Equilibrium (LTE), in which the ^^CO and C^^O lines are assumed to reach the 
same excitation temperature as ^^CO line. The ^^CO (1-0) line is optically thick, so we 
can derive the kinetic temperature directly from its brightness temperature. The kinetic 
temperature is found to be 12 - 14 K as listed in Table [5] and agrees with that derived from 
ammonia. 

The ^^CO and C^^O line intensity ratio indicates the different abundances and optical 
depths in the cloud. As shown in Fig. [31 the ^'^CO to C^^O ratios in regions with strong 
emission are low (about 1.4-2.5), while the ratios of fain t emission areas are high and close to 



the lo cal interstellar p3c0]/[CiS0] ratio of about 7.3 (IWilson fc Roodlll994] : iTevssier et al. 



2OO2I ). The observation of v ariable line ratios as a function of intensities is consistent with 
that of iTeyssier et al. J2002h . 
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3.3. Line widths and profiles 

The NH3 (1,1) hne widths in the cloud vary f rom 0.4 to 2.8 km s~^, which is consisten t 



with the value reported from massive dense cores (iBensen &: Myerslll989t iPillai et al.ll2006l ). 
The distribution of ammonia line widths is shown in Fig. HJ The line widths in the peripheral 
regions (>1.3 km s~^) are generally larger than those in the clump centers (0.7 - 1.4 km s^^). 
The line width of Core 6 is influenced and broadened by the extended part of core 1 and 
double peaked spectra appear in the central part of core 6 (shown in Fig. [5l P6). For cores 
2, 4, 5 and 6, the (2,2) line widths are larger than the (1,1) line widths, which s uggests that 



the ( 2,2) line may not trace the same gas as the (1,1) line toward these cores (jPillai et al. 



20061). 



As seen in Fig. |5l the CO lines show multiple velocity components. Most of the ^^CO 
spectra in the cores are blended with a component around 5.5 km s~^, which also appears in 
some C^^O spectra. This compon ent is consistent w i th th e extended part of another cloud 
to the west of core 1, reported by iFeldt &: Wendkerl (Il993l ). A few of the ^^CO spectra are 



fiat topped, indicating that the line is optically thick at these locations. As with the ^^CO 
line profiles, the C^^O spectra are also slightly asymmetric, especially in core 2. The line 
widths over the core regions range between 2.3 and 4.6 km s~^ for ^'^CO and between 1.5 
and 2.7 km s'^ for C^^O. 

The thermal line width, [(8 ln2)A;Tfcm/m]^''^ {k is the Boltzmann constant and m is the 
mean molecular mass), at a kinetic temperature, Tkin, of 12 K is 0.18 km s~^ for NH3 and 
0.14 km s~^ for ^'^CO and C^^O. Therefore, the non-thermal line widths are significantly 
larger than the thermal line widths and near the observed line widths. This suggests that 
non-thermal broadening mechanisms (rotation, turbulence, etc) play a dominant role in 
producing the observed line profiles. 

The HCO+ (1-0) spectra as shown in Fig. [6] display three types of profiles: an asym- 
metric, double-peaked shape in the northern C^^O condensation (e.g., B, C, PI, P2, etc.), 
spectra with a blue shoulder in the far north (e.g.. A, P6, etc.) and a single blue-shifted peak 
in the south (e.g., D, etc.). The single peaked C^^O line appears to bisect the HCO"*" profiles 
in velocity, indicating that the HCO"^ lines are likely optically thick and self-absorbed. 

The blue-shifted dip due to self-absorption may indicates that this occurs in an expand- 
ing envelope. By fitting Gaussian line profiles to these spectra and comparing th e peak veloc- 



ity of HCO"'" and C^^O spectra, one can roughly estimate an expansion velocity (lAguti et al. 



2003). We measure a velocity of 6V = \vp^^^{C^^O) - t;peak(HCO+)| ^ 1.05 km s-\ This 
value is less than the C^^O line width at corresponding positions but is supersonic. On the 
contrary, the blue peaked profile in the south can be associated with infall asymmetry and 



- 8 - 



may suggest infall motions of the outer cloud material. However, due to the low signal to 
noise ratio of the HCO"'" spectra nearby and the relatively large beam sizes involved, this 
motion can not be confirmed at all positions with blue-shifted peak. In order to better un- 
derstand the spatial distribution of the HCO"*" line profile, we have compared the velocity of 
peak emission in the HCO"'" and C^^O lines for each location in the northern condensation, 
and mapped the velocity difference, 6V = fpeak(C^^O) — WpeaklHCO"*") (Fig. [6]). It can be 
seen that the outflow asymmetry (coded blue) dominates the central and northern regions 
of the cloud, while the infall asymmetry (coded red) appears on the southern edge. 



3.4. Velocity structure 

From Table [3l we note that the line width of core 2 is rather small (0.95 km s^^) but 
remarkably increases to 1.2 km s~^ when the spectra are averaged over the core area, which 
reflects the existence of velocity gradients in the core. This is also seen in the NH3 (1,1) 
channel maps shown in Fig. [7] for velocities from —1.0 to 3.0 km s~^ where different clumps 
appear at separate velocities. The clump structure is complex and covers a wide velocity 
range. In the lower velocity channels (—1.0 to 0.0 km s~^), the gas is concentrated in a clump 
in the east of northern condensation, while the higher velocity channels (0.5 to 2.0 km s^^) 
indicate the clump in the two peak position and a southern clump. 

The CO presents multiple peaks along the line of sight (cf. CO spectra in Fig. |5]). 
The channel maps of ^^CO (1-0) and C^^O (1-0) (Fig. [8]) show that the different velocity 
components present various shapes and orientations. C^^O shows a velocity structure that is 
similar to that of ammonia, shown in Fig. [71 The velocity increases from the position of core 
1 at about 1 km s^^ to 2 km s~^ towards the north, which can be seen in position-velocity 
maps (Fig. |9]) and the ^^CO, C^^O and NH3 channel maps. 

To estima te the magnitude and d irection of the velocity gradient in each core, we adopt 



the method of iGoodman et al.l (Il993[ ) by fitting a linear velocity gradient to the projected 



velocity field along the line of sight. Thus the observed velocity wlsr can be expressed as 

^LSR = vq + ^^Aa sin 9<y -|- $f A5 cos 0c^, 

where ^ is the magnitude of the velocity gradient, Aa and A6 represent offsets in right 
ascension and dechnation in arcseconds, Vq is the systemic velocity of the cloud, and 9^ is 
the direction of increasing velocity, measured east of north. We carry out a least-squares fit 
to the two dimensional velocity field of ^^CO (1-0), C^^O (1-0), NH3 (1,1) and NH3 (2,2) 
weighted by I/ctlsr, where o"lsr is the uncertainty in fitting fLSR- In Table |6l we list our 
fitting results and O^^) with their formal errors (cr<y and ctq^). We exclude fits that fail 
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the 3(7 criterion > Scr^y mentioned by ( iGoodman et al.lll993l )) or those with a random 
velocity field as seen from the velocity map. 

Both cores 1 and 6 present significant velocity gradients in CO and NH3 lines. The 
magnitude of the velocity gradient is greater in NH3, and the direction of the gradient is 
somewhat different from that seen in CO. In cores 3, 4 and 5, the velocity gradient is detected 
only in NH3 (1,1). The NH3 gas associated with core 2 also exhibits a clear velocity gradient 
in Fig. m The peak lies on a velocity ridge with a gradient of 3.15 km s~^ pc~^ to the east 
and 1.93 km s~^ pc~^ to the west. 

The distinct velocity gra dients given by differe nt molecules are probably due to the 
different densities they trace ( IGoodman et al.lll993[ ). Velocity gradients seen in the lower 
density tracers, ^^CO or C^^O, probably arise in the envelope of the northern condensations. 
This gradient is dominated by unresolved clump-clump motions within the envelope, and 
may also be contaminated by the higher velocity component. The higher density tracers 
(e.g., NH3), on the other hand, reveal the gradients within the small clumps themselves. 



3.5. Density and abundance 



The column densities of CO are estimated using the equation given by IScoville et al. 



(]1986f ): 



N 



3k' 



exp 



hpj\ T^^ + hu/6k{J +1] 



exp(— /iz//A;Tex) 1 — exp(— r) 



where J is the rotational quantum number of the lower state in the observed transition, /x is 
the permanent dipole moment and r is the optical depth derived from the observed ratios of 
^^CO and ^^CO (or C^^O) emission. The derived physical properties are listed in Tabled We 
also derive the column densities of NH3 for the cores, which range from 15 — 34 x 10^^ cm~^ as 
shown in Table O Since ^^CO is not optically thin in some cores, we estimate the H2 column 
densities based on the C^^O column densities wher e a factor iV(H2)/A^(C^^0) = 7 x 10^ was 
adopted (IFrerking et al.l Il982l : iKramer et al.l Il999l ) . This gives H2 column densities around 



~ 10^^ cm ^ as listed in column (2) of Table [71 

With the H2 column densities derived from CO emission, we determine the ammonia 
abundances and list them in Table [71 The ammonia abund ance in the s ix pe ak positions is 
3 — 5.5 X 10~^, which agrees with the parameters reported by lPillai et al.l (120061 ). though their 
abundances are derived from the dust emission of SCUBA. When averaged over the cores, 
the abundances are reduced to 2.3 x 10"*^. The ammonia abundance and the correlation of its 
spatial distribution with the morphology of the dust emission is consistent with a chemical 
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model given by iBergin &: Langerl (119971 ). in which NH3 is not depleted. We then compare 
its abundance with that of ^^CO toward cores in Fig. [TOl They differ slightly from each 
other with core 1 showing the highest abundance among the cores. This reveals the complex 
chemical environment in the cloud. 

We obtain the gas mass of cores by integrating the total H2 column density over the 
cores and show the results in Table [HI The mass range from 60 Mq to 250 Mq. The gas 
density can be derived by assuming the cores to be spherical and dividing the mass by the 
volume of the cores. An average density of ~ 1.1 x 10^ cm~'^ was derived as shown in Table El 

The total H2 column density can also be related to the 1.1 mm flux by assuming the 
continuum emission to be optically thin. After convolving the data to the resolution of the 
100 m Effelsberg (40") and 13.7 m Delingha (60") beams, the dust based column density is 
given by 

mm 



2)1.1 mm ^ j3 frp \ 1 

"bcam-D!/UdJ'tl.l ramTTl 

where f^beam is the telescope beam solid angle, B^iTd) is the Planck function at the dust 
temperature (assumed to be equal to the gas temperature deduced from NH3), i mm 
is the dust opacity at 1.1 mm, and m is the mean molecular niass. H ere we adopt ki.i mm 
to be 0.0114 cm^ g~^ (jOssenkopf fc Henninglll994l : Enoch et al.ll2008l ) and assume a gas to 
dust mass ratio of 100. The results have the same magnitude as those derived from C^^O 
and are listed in Table [71 We also compare our molecular emission with the BGPS survey, 
as shown in Fig. [TTl The integrated intensity of NH3 (1,1) shows a better correlation with 
the 1.1 mm dust emission compared to C^^O as mentioned previously. 

The diversity in column densities derived from molecular gas and dust may be due to 
different structures and components they trace. The CO column density, because of the 
low density it traces, is highly affected by the structure of the cloud along the line of sight, 
especially at the positions where clumps with similar velocities blend together. Also the 
molecular observations are affected by high opacity or chemical depletion at high densities. 
On the other hand, calculations related to the dust emission are based on the value of dust 
opacity and assumptions about the dust to gas ratio. The d ust opacity changes betwee n 
different dust models and in environments of different density (jOssenkopf fc Henninglll994l ). 
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DISCUSSION 



4.1. Evidence of massive star formation 



In this section, we will compare our results with those from other low-mass star forming 
regions, and discuss the evidence for massive star formation and the evolution ary state of 
t he cl oud. The mean intrinsic line width for NH3 (1,1) in the samples given b vjjiiina et al. 
19991) (most of the objects are low-mass cores) is 0.5 km s~^. A recent work by lFoster et al. 



(I2OO9I ) also reported a sm all observed line wid th of 0.3 km s ^ for low mass NH3 cores in the 



Perseus Molecular cloud. ICrapsi et al.l ( 120051 ) undertook a survey toward low-mass starless 
cores and also reported a narrow line width of 0.5 km s~^ for C^^O (1- 0) . A broader lin e widt h 
of 1.2 km s~^ for C^^O (1-0) is derived from a JCMT observation by |j0rgensen et al.l (120021 ). 
which is still less than those given in our observed region. Meanwhile, the dominant non- 
thermal line widths in our cloud indicate more non-thermal support than those in low-mass 
star forming regions. Although line widths can be broadened by active outflows in low-mass 
protostellar regions to 2 km s~^ for C^^O (1-0) JSwift fc Welch! boosl ) and 1.2 km s'^ for 
NH3 (1,1 ) (IRudolph et al.ll200ll ). at this stage, an active o utflow has not de veloped yet. 
The NH3 line widths here are comparable to those IRDCs of iPillai et al.l ( 120061 ). Moreover, 



our derived column densities of NH3, a few times 10^^ cm ^, a re higher than those of the 



low-mass cores (~ 10^^ — 10^^ cm ^) given by [ Suzuki et al.l (119921 ). and are also comparable to 



those of the massive IRDCs given by lPillai et al.l (120061). In addition, the supersonic motions 
of the envelope and the large velocity gradients in our cloud also indicate a more dynamic 
motion compared to low-mass cores. Therefore, we suggest that MSXDC G084. 8 1—01.09 is 
potentially at an early evolutionary state of massive star formation. 



4.2. Gravitational stability of the cores 

The expanding envelope detected from the line profiles described in §3.31 lead us to 
consider the gravitational stability of MSXDC G084. 81— 01.09. By using the molecular mass, 
M Molecular, derived from C^^O (1-0) and the core radius R determined from the NH3 (1,1) 
maps, we obtain an escape velocity {^2GM/R) which ranges from ~ 1.9 km s in core 6 to 
a maximum of ~ 3.2 km s~^ in core 2 (Table |8|). The three-dimensional velocity dispersion 
of the cores provided by the NH3 line width lies in the range of 1.2 - 2.4 km s~^ which 
is less than the escape velocity in each core. However, the C^^O velocity dispersions (2.5 
- 4 km s~^) exceed the escape velocity in some regions of the cloud, especially for core 6, 
which has a small escape velocity. This may lead to instability for core 6. Additionally, an 
expansion speed of 1.0 km s~^ in the envelope of the core (cf. §3.3p . though significantly less 
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than the escape velocity, may aggravate the situation. 

To examine the gravitational stability of the cores, it is useful to calculate the virial 
mass of the cores as well as the virial par ameter, the ratio between the virial mass and 



core mass. The virial parameter given by iBertoldi fc McKed ( 1l992l ) can be expressed as 
a = 5a'^R/GM, where R is the radius of the core, a = A/3/(81n2) x FWHM, and M 
is the core mass. The virial mass and virial parameter can be found in Table El The 
average virial parameter is about 1.1 which suggests that most of the cores are virialised. 
However, the virial parameters of cores 2 and 6 deviate significantly from unity. This can be 
caused by a number of factors. One explanation is the potential for the presence of spatially 
unresolved stellar clusters in the cores, which can be identified from the color temperature 
of the dust continuum map. For Myiriai > ^Molecular, beam filling factors smaller than unity 
or streaming motions may adversely affect the estimate and may lead to an underestimated 
column density or an overestimated line width. Alternatively, it is po ssible that the cores 
may not actually be in virial equilibrium, and may be transient entities (IBallesteros-Paredes 



20061 ). which is probably the case in core 6. 



We derive a Jeans mass greater than 170 Mq for each core using the equation Mj^ans = 
17.3 Tfcj„^'^n~°'^MQ, where Tfcj„ is the kinetic temperature derived from NH3 data, and n 
is the average density as listed in Table. [H We note that the masses of cores 1 and 2 
are comparable to their Jeans masses, but the other cores are significantly less massive. 
However, considering the clumping in cores 1 and 2, it is still possible for them to form 
several individual protostars. Therefore, we suggest that cores 1-5 are gravitationally 
bound, among which cores 1 and 2 are marginally stable against collapse. However, core 6 
is probably transient. 



4.3. Comparison with more evolved clouds 



Our dark cloud has a lower kinetic temperature than other evolved clouds. IWu et al. 



( 120061 ) reported a mean temperature of 19 K in massive water maser sources ex c luding 



known HH regions, and in ultra-compact HII (UC HII) regions. IChurchwell et al.l (119901 ) 
give a significantly higher temperature spread from 15 K to >60 K for UC HH regions. The 
temperature of the cloud is a good tracer to determine its evolutionary stage in the context 
of massive star formation. The low temperatures in our cloud suggests that it is not an 
active high- mass star forming region yet. 

The C^^O line widths are typically a factor of two to five times smaller than the line 
widths in more evolved massive star forming regions traced by other molecules with similar 
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densities. The NH 3 li ne widths we find are comparable to the massive cores reported by 
Pillai et al.l ( 120061 ) and IWu et al.l ( 120061 ). but are smaller than those associated with water 
masers or are in proximity to UC HII regions. More specifically, the line widths of our cloud 
are similar to those of the most quiescent NH.c; co r es wi th no methanol maser or 24 GHz 
continuum emission described by iLongmore et al.l (120071 ) and smaller than those in more 



evolved cores with maser or continuum emission. Our cloud also presents similar NH3 line 
width to those bright-rimmed clouds which are undergoing r ecently initiated sta r formation 
and being subjected to intense levels of ionizing radiation ( iMorgan et al.ll2010l ). After re- 
moving the thermal broadening contributions to the line widths, the velocity dispersions 
associated with turbulence are supersonic and lie between those of the sources triggered by 
the radiatively driven implosion and the non-t riggered ones. For the more evolved clouds 
in the UC HII phase, IChurchwell et al.l (Il990[ ) reported significantly larger line widths of 
~ 3 km s~^. The large line widths in more evolved stages are likely due to a combination of 
warmer temperatures and broadening from dynamics such as outfiows. 

The HCO+ (1~0) spectra associated with hyp er-compact H II r e gions have a consider- 
ably larger line width (15-60 km s~^) reported bv IChurchwell et al.l (120101) . For our cloud, 



our H CO^ line profiles are similar to the samples of iBrand et al.l ( 1200 ll ) and ICesaroni et al. 



(Il999h. which are all in the pre-UC HII stage of star formation, though some sourc es in 
Cesaroni et al.l (Il999[ ) are proven to be more evolved than those in iBrand et al.l ( I2OOII ). We 
find similar line widths and profiles as theirs, with a self-absorption dip between the blue 
and red peak. However, our cloud is colder than the pre-UC HII sources and do not have 
any associated IRAS point source. 



The optical depths of N H3 (1,1) are comparable to the cores given by lPillai et al.l (|2006[ ) 
and iLongmore et al.l ( 120071 ). Our core 4, 5 and 6 have a lower optical depth than the NH3 
sources (~ 2.7 — 2.9) associated with methanol masers or 24-GHz continuum emission in 
Longmore et al.l ( 120071 ) . while core 1, 2 and 3 near the cloud center gives higher optical 
depth but still lower than those samples (~ 3.1) with only NH3 association. However, due 
to the large uncertainty, without further observation, it would be premature to classify these 
sources into more detailed evolutionary stages based on optical depth. Our cloud density 
is about an order of magnitude lower than the typical val ue in more evo l ved rn assive star- 
forming regions: an average density of 10^ cm~^ is given by lBeuther et al.l (|2002[ ) using LVG 
calculations in clo uds prior to building up an UC HII region, and a similar value is given by 
Pillai et al.l (120071 ) in clouds harboring a UC HII region. 



Our NH3 (1,1) core sizes lie between themean value of 0.28 pc given by ILongmore et al. 



( I2OO7I ) and 0.57 pc given b y iPillai et a. 



(|2006[) . but are significantl y smaller than th e mean 



size of 1.6 pc reported by IWu et al.l (120061 ). The small values of ILongmore et al.l (|2007| ) 
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probably result from their small beam size (11")) but the large values of IWu et al.l (|2006[ ) 
are likely a result of evolution. All these aspects suggest that our cloud is in an early 
evolutionary stage of massive star formation. These quiescent cores are likely the candidates 
for pre-protostellar cores. 



4.4. Rotation in the cores 



One of t he possible generators o f the velocity gradients mentioned in §3.41 is rotation 
in the cores. iGoodman et al.l (Il993l ) has discussed the relation between rotation and the 
geometric properties of the cores. They found no causal relation between velocity gradient 
direction and core elongation, nor any relationship between the magnitude of the gradient 
and core shape, on the size scale of 10^^ cm. An examination of the gradients in our clouds 
also reveals no obvious relation of this kind. In view of the early evolutionary state and 
the dominant role of non-thermal line broadening in the cloud, we believe that the velocity 
gradients in the cloud are more affected by stochastic processes, such as fragmentation, 
collisions, and nonuniform magnetic fields rather than ordered motions such as rotation. 



We use the parameter (3 as defined by lGoodman et al.l (119931 ) to compare the rotational 
kinetic energy to the gravitational energy. Thus /3 can be written as 

{l/2)Iio^ _ Ipco^R^ 
^ ~ qGMyn ~ 2q GM ' 

where I is the moment of inertia given by / = pMB?, qGM'^/R represents the gravitational 
potential energy of the mass M within a radius R, and w = ^/ sini, where i is the inclination 
of the core along the line of sight. We assume p/q = 0.22 as for a sphere with an density 
profile and sini = 1. Our f3 values as listed in Table [6] are consistent with the result of 



Goodman et al.l ( 1l993l ) that most clouds have (3 < 0.05. The small values of /3 show that 
the effect of rotation is not significant in maintaining the overall dynamical stability for the 
cloud. The results also indicate that these cores are unlikely to experience instabilities driven 
by rotation (e.g., bars, fission, or rings), if the magnetic fields are not taken into account. 



5. Summary 

The multiple molecular line observations of the infrared dark cloud MSXDC G084. 81— 01.09 
were analyzed. The results are summarized as below. 

The cloud has a low gas temperature of 12 K, a column density of 10^^ cm~^ and 
masses of the cores range from 60 to 250 Mq. The spatial distribution of the ammonia 
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emission correlates well with the mid-infrared absorption and millimeter dust emission. All 
6 ammonia cores are associated with their corresponding dust peaks. The line width of 
1 km s~^ for ammonia indicates the dominant role of non-thermal broadening in the cloud. 
The abundances of ammonia range from 2 — 3.5 x 10^^. These facts show that the cloud is 
a site of potential massive star formation at a very early evolutionary stage. 

""^^CO and C^^O trace a more extended cloud structure. Together with the HCO"*" emis- 
sion, we detect an expanding envelope in part of the cloud with a expansion velocity of 
1 km s~^. Five of the cores are gravitationally bound (four are virialised) with Jeans masses 
comparable or exceeding their molecular mass except a possible unbound transient core. We 
found velocity gradients in the cores using different molecular tracers. The different direc- 
tions and magnitudes of the gradients using the different tracers indicate clumping motions 
on different scales. 

The observations presented here are part of a larger program designed to search and 
investigate potential massive star forming regions. Higher resolution observations with inter- 
ferometers would be needed to study turbulence and fragmentation inside our cores. Addi- 
tional submillimeter mapping observations could derive the dust temperature and structure 
associated with the gas. We could also use radio polarization observation to study the re- 
lationship between magnetic field and cloud structure at the early phrases of massive star 
forming. 

This work is made based on observations with the 100-m telescope of the MPIfR (Max- 
Planck-lnstitut fiir Radioastronomie) at Effelsberg and Dclingha 13.7-m telescope of Purple 
Mountain Observatory. This work is also based on observations made with the Spitzcr 
Space Telescope, which is operated by the Jet Propulsion Laboratory, California Institute 
of Technology under a contract with NASA. The authors appreciate all the staff members 
of the observatories for their help during the observations. This work was supported by 
the Chinese NSF through grants NSF 11073054, NSF 10733030, NSF 10703010 and NSF 
10621303, and NBRPC (973 Program) under grant 2007CB815403. 
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Fig. 1. — Spitzer MIPS 24 fim image overlaid with an NH3 (1,1) contour map of integrated 
intensity from the MSX dark cloud G084. 81— 01.09. The spectra were integrated over —2.5 
to 4.5 kms~^. The peak contour flux is 8 Kkms^^ with contour interval 1 Kkms~^. The red 
star indicates the ioii i zing s tar of W80 complex, 2MASS J205551. 25+435224. 6, reported by 
Comeron fc Pasaualil tood) The cross indicates the MYSO candidate, G084.7847-01.1709, 
identified bvlUrgiiiart et al. (2009). The green triangles indicate the star clusters detected by 
Camberesy et al.l (120021 ). The white circle shows the beamsize of the NH3 (1,1) observation. 
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Fig. 2. — CSO BOLOCAM 1.1 mm continuum image overlaid with an integrated intensity 
contour map of observed molecular lines, (a) NH3 (1,1) integrated over —2.5 to 4.5 kms~^, 
with a bold line at 50% of maximum intensity of 8.7 Kkms~^ in steps of 9%. (b) ^^C0(1— 0) 
integrated over —0.5 to 2.5 kms~^, with a bold line at 70% of maximum intensity of 38.2 
Kkms~^ in steps of 5%. (c) 0^*^0(1-0) integrated over —0.5 to 2.5 kms~^, with bold 
line at 70% of the maximum intensity of 8.3 Kkms~^ in steps of 6%. (d) HC0^(1— 0) 
integrated over —2.5 to 4.5 kms~^, with bold line at 70% of the maximum intensity of 7.7 
Kkms~^ in steps of 8%. Four peaks in the HCO+ map are design ated as A-D. The cros s 
indicates the MYSO candidate, G084.7847-01.1709, identified bv lUrquhart et al.l feoooh . 
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Fig. 3.— Tmb of i3CO(l-0) vs. that of 0^^0(1-0) for G084.81-01.09. The data are 
extracted within a velocity range of —1.5 to 3.7 km and smoothed to 0.2 km s~^ spectral 
resolution. Each dot represents the Tmb of the same channel at same position for ^^CO 
and C^^O. Data points from the northern and southern clumps are shown in black and grey 
colors respectively. The dashed-linc represents the local ISM ratio. Note that the northern 
and southern clumps have a different trend for the intensity ratio at high intensities. 



- 23 - 




20''57'"15.0^ SV^O.O^ 56'"45.0^ 56'"30.0^ 20''57'"15.0^ 57'"0.0^ 56'"45.0^ 56'"30.0^ 

Right Ascension (J2000) Right Ascension (J2000) 



Fig. 4. — Left panel: Color-coded map of the NH3 (1,1) full width at half maximum (FWHM) 
line width in the MSX dark cloud G084. 81— 01.09 overlaid on a contour map of the integrated 
intensity. Right panel: Color-coded map of the local standard of rest velocity (Visr) overlaid 
on a contour map of the integrated intensity. Only spectra with signal greater than 5a 
(FWHM map) and 3cr (Visr map) are shown. The circle in the bottom right shows the 
beamsize. 
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Fig. 5. — Spectra of G084. 81— 01.09 toward the six peak positions (labeled as P1-P6). The 
first two rows show NH3 (1,1) and (2,2) spectra. The NH3 (2,2) lines are moved upward in 
each panel for clarity. The main (1,1) line feature of P6 was extracted and drawn beside the 
spectra. The next two rows show the ^^CO (1—0) and C^^O (1—0) spectra. In all cases, the 
stronger line belongs to ^^CO. The peak positions are referred to in Table |2l 
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Fig. 6. — Left: HCO"'" (solid) and C^^O (dashed) spectra at the positions of the 4 NH3 cores 
and the 4 HCO"*" peaks (indicated by red characters in Fig. [2]) from north to south. The 
vertical dashed lines indicate the Visr of the NH3 (1,1) line, obtained from a Gaussian fit. 
Right: Map of SV, the velocity difference between the peak of C^^O and HCO"'" overlaid on 
a C^®0 contour map of the region. 
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Fig. 7. — NH3 (1,1) channel maps for G084. 81— 01.09. The central velocity of each channel, 
in kms~^, is marked on the top- left corner for each map. The reference position is at 
20''57™00.9^+43°44'15.4". 
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Fig. 8.— ^3CO(1-0) channel maps (upper panels) for G084.81-01.09 and 0^*^0(1-0) chan- 
nel maps (lower panels) for G084. 81— 01.09. The central velocity of each channel, in kms~^, 
is marked on the top-left corner for each map, and the plus signs indicate the positions of 
6 cores identified from NH3 (1,1) map (see Fig. [T]). The reference position is the same as in 
Fig.El 
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Fig. 9. — The position- velocity map of ^^CO (top), C^^O (second row), HCO^ (third row) 
and NH3 (1,1) main component (bottom) along the Declination axis for the MSX dark cloud 
G084.81-01.09, at Right Ascension of 20^^56"^47.1" (left panel with offset A(5=-150") and 
Right Ascension of 20^^57^03. 7*^ (right panel with offset A5=30"). The declination reference 
coordinate is the same as in Fig. [3 
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Fig. 10. — Abundances relative to ^'^CO toward peak positions (Table [2]). The relative 
abundances were shown in different grey scale. 
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Fig. 11. — Relation between the NH3 (1,1) and C^^O (1-0) integrated intensity and the 
1.1 mm flux density. Continuum data were smoothed to corresponding resolution. The 
uncertanties of integrated i ntensity are the rms nois e of the spectra and the error in dust 
flux is calculated following iRosolowsky et al.l ( 120091 ). The solid line represents the least- 
square fit of the data weighted in both coordinates. 
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Table 1. Observation parameters 



Line 


2^0 


Telescope 


During 


HPBW 


6u 


5v 






(GHz) 






(") 


(KHz) 


(km s^-') 


(K) 


NH3 {J,K=1, 1) 


23.694496 


Effclsbcrg 100m 


2007 Dec. 


40 


19.5 


0.24 


30-150 








2008 Feb. 


40 


30.5 


0.38 


30-45 


NH3 (J, K=2, 2) 


23.722633 


Effelsberg 100m 


2007 Dec. 


40 


19.5 


0.24 


30-140 








2008 Feb. 


40 


30.5 


0.38 


30-45 


NH3 (J, K=2,, 3) 


23.870130 


Eflfelsberg 100m 


2007 Dec. 


40 


19.5 


0.24 


30-140 








2008 Feb. 


40 


30.5 


0.38 


30-45 


NH3 (J, K=i, 4) 


24.139417 


EfEclsbcrg 100m 


2007 Dec. 


40 


19.5 


0.24 


30-160 


i2CO (J=l-0) 


115.271204 


DLH 13.7m 


2008 Mar.-Apr. 


60 


142.1 


0.37 


200-400 


"CO (J=l-0) 


110.201353 


DLH 13.7m 


2008 Mar.-Apr. 


64 


41.7 


0.11 


200-400 


cisq (J=1-0) 


109.782183 


DLH 13.7m 


2008 Mar.-Apr. 


64 


42.1 


0.11 


200-400 


HCO+ (J=l-0) 


89.188530 


DLH 13.7m 


2008 May.-Jun. 


~78 


42.1 


0.14 


270-420 



Note. — vq is the rest frequency of the hne, HPBW is the half power beam width of the telescope, 5v and 5v 
represent the frequency and velocity resolutions 
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Table 2. Cores detected in the NH3 observations of G084.81-01.09 



Peak 


a 


5 


Transition 


14sR 








FWHM 


''"main 






(J2000) 


(J2000) 




(kms-l) 




(K) 




(kms-i) 




(Kkms-l) 


1 


20 56 47.1 


+43 44 45 


1-1 


1.25+0.01 


3, 


.84+0, 


,24 


1.38+0.03 


3.63+0.20 


8.73+0.20 








2-2 


1.30+0.04 


1, 


.14+0. 


,19 


0.96+0.14 




1.39+0.19 


2 


20 56 47.1 


+43 43 15 


1-1 


1.23+0.01 


4, 


.98+0, 


,37 


0.75+0.02 


3.61+0.24 


7.26+0.19 








2-2 


1.20+0.03 


1, 


,51+0, 


,17 


0.88+0.09 




1.61+0.19 


3 


20 57 03.7 


+43 40 45 


1-1 


1.22+0.01 


3, 


,67+0, 


,36 


0.89+0.04 


3.25+0.29 


6.03+0.19 








2-2 


1.20+0.01 


1 


,07+0, 


,19 


0.63+0.08 




1.09+0.19 


4 


20 57 00.9 


+43 38 45 


1-1 


0.99+0.01 


3, 


,71+0, 


,31 


1.00+0.03 


2.22+0.20 


5.91+0.18 








2-2 


0.84+0.10 


0, 


,58+0, 


,14 


1.40+0.28 




0.78+0.18 


5 


20 57 03.7 


+43 37 15 


1-1 


1.20+0.01 


3, 


,52+0, 


,25 


0.97+0.03 


2.32+0.17 


5.53+0.10'* 








2-2 


1.21+0.06 


0, 


,72+0, 


,09 


1.32+0.21 




1.02+0.11 


6 


20 56 55.4 


+43 47 15 


1-1 


2.01+0.02 


2, 


,51+0, 


,23 


1.40+0.04 


2.08+0.21 


5.02+0.18 








2-2 


2.07+0.14 


0, 


,53+0, 


,12 


1.98+0.34 




1.39+0.19 



Note. — Columns arc peak number, offset position from map center, NIfs (J,K) transition, the local standard rest 
velocity, main beam brightness temperature, full width at half maximum, optical depth and integrated intensity of 
the main line. Values in column (4), (5), (6) and (7) are the HFS or GAUSS fitting results and errors estimated in 
CLASS. Integrated intensities are calculated from —2.5 to 4.5 kms~^ with errors derived from rms. 

^PeaJj 5 was observed in FFT mode with lower noise compared to Peak 1-4 and 6 in frequency switch mode. 
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Table 3. Physical properties of the NH3 cores 



Number Tex Trot"" Tkin^ Ar(NH3)<= AiJi^ Aj^a'^ 

(K) (K) (K) (10"cm-2) (kms-i) (kms-i) 

1 6.7±0.5 11.4±0.6 12.2±0.7 33.5±2.9 1.30±0.01 1.56±0.01 

2 7.8±0.8 11.5±0.5 12.3±0.6 21.1±2.2 0.95±0.01 1.21±0.01 

3 6.5±0.8 11.7±0.7 12.5±0.9 18.3±2.5 0.87±0.01 0.98±0.02 

4 6.9±0.8 10.5±0.7 11.1±0.9 17.5±2.2 1.06±0.01 1.19±0.02 

5 6.6±0.6 11.2±0.5 11.9±0.6 15.3±1.6 l.OliO.Ol 1.05±0.02 

6 5.6±0.8 11.5±0.8 12.4±1.0 16.1±2.2 1.21±0.02 1.49±0.03 



Note. — Columns are peak number, excitation temperature, NH3 rotational tempera- 
ture, kinetic temperature, NH3 column density derived from the optical depth, the mean 
line width and the line width of the mean spectrum of the regions. 

^NHa rotational temperature given bv lHo fc Towne^ l|l9S3l ) as 



-41.5^ In 



-0.282 

Tm(l,l) 



■ In 



r*(2,2,m) , 

1 — L X (1 - 

r*(l,l,m) 



^NHs kinetic temperature Ti^jn 



at ln[l+l.lcxp{-16/Trot)l 



'^NHa column den sity derived from its optical depth via the relation given by 
iBachiller et al. 1 lll987l) as 



Ar(NH3) = 2.784 x 10"rJ(Tcx)Ay X Q/Qi, 

where Q is the partition function and Q/Qi = ie23-4/Trot + 1 + 5^-41. 5/T„t _|_ 
Wg-101.5/T„t _|_ ,,,, 

"^The mean line width is the fitting error- weighted average line width over the core 
regions. 

"^The line width of the mean spectrum is the line width of mean spectrum of the 
spectra inside the core regions. 
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Table 4. CO and HCO"'" line parameters at the positions of the NH3 cores. 



Peak 


Molecule 


Vlsr 




AV 




(7 






(kms-i) 


(K) 


(kms-l) 


(Kkms-i) 


(K) 


1 


12CO (1-0) 


2.31 


10.76 


6.38 


89.21±0.44 


0.20 




"CO (1-0) 


0.96 


8.00 


3.19 


31.44±0.22 


0.22 




C^^O (1-0) 


0.83 


3.04 


2.17 


6.41±0.10 


0.17 




HCO+ (1-0) 


1.44 


1.00 


4.76 


4.99±0.25 


0.21 


2 


12C0 (1-0) 


2.67 


11.02 


6.47 


89.02±0.72 


0.33 




13CO (1-0) 


1.03 


7.48 


3.08 


30.84±0.22 


0.22 




Ci»0 (1-0) 


1.04 


2.75 


1.83 


5.21±0.13 


0.21 




HCO+ (1-0) 


1.41 


1.61 


3.11 


5.83±0.20 


0.17 


3 


12CO (1-0) 


2.21 


9.86 


7.24 


95.16±0.51 


0.23 




13CO (1-0) 


0.90 


5.87 


4.05 


28.53±0.25 


0.25 




Ci»0 (1-0) 


1.08 


2.51 


1.72 


4.49±0.15 


0.24 




HCO+ (1-0) 


-0.52 


1.08 


2.11 


2.44±0.24 


0.20 


4 


(1-0) 


2.64 


9.79 


6.69 


85.58±0.44 


0.20 




"CO (1-0) 


1.11 


5.71 


3.30 


25.15±0.21 


0.20 




Ci^O (1-0) 


1.04 


2.83 


1.62 


4.80±0.10 


0.15 


5 


i^CO (1-0) 


2.80 


9.76 


6.34 


77.12±0.43 


0.20 




"CO (1-0) 


1.32 


5.57 


3.14 


23.28±0.21 


0.21 




C"0 (1-0) 


1.24 


2.47 


1.50 


3.99±0.11 


0.17 


6 


12CO (1-0) 


2.47 


10.89 


5.85 


77.47±0.45 


0.20 




"CO (1-0) 


1.67 


8.34 


2.70 


28.82±0.23 


0.23 




C"0 (1-0) 


1.57 


2.62 


2.08 


5.18±0.10 


0.16 




HCO+ (1-0) 


2.08 


1.35 


3.33 


4.91±0.19 


0.16 



Note. — l^Qo, l^CO, C^^O and HCO+ line parameters are calculated within 
the velocity range from -3.5 to 6.8 kms^^, -2.7 to 7.7 kms~^, -1.0 to 3.0km 
and 0.0 to 9.0 km s^^, respectively. The HCO+ line was not observed towards 
cores P4 and P5 due to low signal-to-noise ratio. 
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Table 5. 


Physical properties derived from CO. 




12CO 


13 


CO 


C 


18o 


Peak 




r 


N 


T 


N 




(K) 




(IQl^cm-^) 




(IQl^cm-^) 


1 


1/1.17±0.20 


1.26±0.09 


58.36±2.12 


0.35±0.03 


8.91±0.31 


2 


14.43±0.33 


1.18±0.09 


56.27±2.32 


0.31±0.03 


6.35±0.33 


3 


13.26±0.22 


0.90±0.07 


45.31±1.55 


0.34±0.03 


6.37±0.32 


4 


12.90±0.19 


0.91±0.06 


39.87±1.18 


0.39±0.03 


6.63±0.30 


5 


12.25±0.20 


1.01±0.08 


37.02±1.29 


0.33±0.03 


4.70±0.26 


6 


14.30±0.21 


1.35±0.10 


54.75±2.19 


0.30±0.02 


7.88±0.29 



Note. — First column is core number. Column densities are calculated under the 
assumption of local thermodynamic equilibrium (LTE), in which ^^CO and C^®0 have 
the same excitation temperature as that of optically thick ^^CO. 
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Table 6. Velocity gradients detected in the cores 



Core 


Molecule 




!y ± (Tc# at 580 pc Beg 


■ ± o"e,^ 


/3 






(kms~^) 


(ms"-'^ arcsec^-'^) (kms"-*^ pc"-"^) (deg 


; E of N) 





1 


i^CO (1-0) 


1.01±0.003 


4.25±0.096 


1.51 


-8.3±0.8 






C^^O (1-0) 


0.91±0.007 


3.69±0.205 


1.31 


-12.1±2.0 






NH3 (1,1) 


1.15±0.014 


10.76±0.339 


3.83 


-37.2±1.2 


1.3E-2 




NHg (2,2) 


1.28±0.024 


8.76±0.927 


3.12 


-50.5±4.4 




3 


NH3 (1,1) 


1.19±0.012 


3.69±0.337 


1.31 


165.0±4.0 


2.3E-3 


4 


• NH3 (1,1) 


l.OliO.Oll 


4.32±0.419 


1.54 


-64.5±4.0 


3.6E-3 


5 


• NH3 (1,1) 


1.12±0.013 


2.31±0.497 


0.82 


-123.8±8.5 


6.2E-4 


6 


"CO (1-0) 


1.66±0.004 


2.87±0.188 


1.02 


16.5±2.2 






CISQ (1-0) 


1.54±0.010 


4.19±0.507 


1.49 


12.6±4.5 






NH3 (1,1) 


l.SlitO.OlS 


9.45±0.772 


3.36 


90.2±3.9 


1.5E-2 




NH3 (2,2) 


1.88±0.048 


14.92±2.772 


5.31 


45.6±7.6 





Note. — Columns are core number, molecules fitted, systemic velocity, magnitude of the velocity 
gradient, at cloud distance, direction of increasing velocity (measured east of north), and parameter 
f) given in i|4.4l Errors quoted are 1 cr uncertainties. 



Table 7. NH3 abundance in the cores 



Peak Af(H2)ci8o XNH3 A^(H2)i.i2mm Xnhs X in core 
(1022cm-2) (10-») (1022cm-2) (IQ-*) (IQ-**) 



1 6.23 5.4 5.36 6.3 2.41 

2 4.45 4.7 5.54 3.8 2.74 

3 4.46 4.1 2.85 6.4 2.46 

4 4.64 3.8 2.16 8.1 2.24 

5 3.29 4.7 0.96 15.9 1.88 

6 5.52 2.9 1.25 12.9 2.04 



Note. — The columns show the core number, H2 column density and 



NH3 abundance x based on C^^O observations (column 2 and 3) and 
1.1 mm dust flux (column 4 and 5), and the average NH3 abundance over 
the cores. 
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Table 8. Masses of the NH3 cores 



Region R n(H2) Mjeans -Wyirial A^Moleoular « ^escape 

(arcsec) (10*cm-3) (Mg) (Mq) (Mq) (kms"!) 

1 69.0 13.8 198 206 208 0.99 3.03 

2 74.9 13.2 205 121 255 0.41 3.22 

3 69.0 8.7 259 92 130 0.71 2.40 

4 60.1 10.4 198 119 104 1.15 2.30 

5 46.7 16.6 174 85 77 1.10 2.24 

6 49.7 11.0 228 129 62 2.08 1.95 



Note. — Columns are core number, size of the cores, H2 densities derived from C^^O, 
Jeans mass, virial masses from NH3 and gas mass, virial parameter a, and escape velocity. 



